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Analytical approximations for (ip 2 ) and (T£) of a quantized scalar field in static spherically sym- 
metric spacetimes are obtained. The field is assumed to be both massive and massless, with an 
arbitrary coupling £ to the scalar curvature, and in a zero temperature vacuum state. The ex- 
pressions for (ip 2 ) and (T^) are divided into low- and high-frequency parts. The contributions of 
the high-frequency modes to these quantities are calculated for an arbitrary quantum state. As 
an example, the low-frequency contributions to (ip 2 ) and (T£) are calculated in asymptotically flat 
spacetimes in a quantum state corresponding to the Minkowski vacuum (Boulware quantum state). 
The limits of the applicability of these approximations are discussed. 
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■ I. INTRODUCTION 

CD '. 

The interest in the vacuum polarization effects in strong gravitational fields is connected mainly with investigations 
Q\ • of the early Universe and the construction of a self-consistent model of black hole evaporation. The main objects to 
calculate from quantum field theory in curved spacetime are the quantities (f 2 ) and (T^) where <p is the quantum 
field and is the stress-energy tensor operator for tp. The latter quantity is of particular interest as a source term 
in the semiclassical Einstein field equations: 



G^ = 8nG(T^). (1) 
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[ Except for very special spacetimes, on which quantum matter fields are propagated, and for boundary conditions with 
■ a high degree of symmetry (see, for example, it is not possible to obtain exact expressions for these quantities. 
(*C) ' Numerical computations of these quantities are as a rule extremely intensive 0, 0- 13 ■ Thus it is useful, when 
possible, to have analytical approximations to (tp 2 ) and (T{f). One of the most widely used techniques to obtain 
information about these quantities is the DeWitt-Schwinger expansion |(J. It may be used to give the expansions for 
(ip 2 ) and (T£) in terms of powers of mL where m is the mass of the quantized field and L is the characteristic scale of 
change of the background gravitational field. For conformally coupled massless fields approximate calculations have 
also been made. For (TJf) in static Einstein spacetimes (i? M „ = Agy,^) these include the approximations of Page, 
' Brown, and Ottewill 0,11, El- These results have been generalized to arbitrary static spacetimes by Zannias [ljj. A 
different approach to the derivation of approximate expressions for ( y>) and (Tjf) for conformally coupled massless 
fields in static spacetimes has been proposed by Frolov and Zel'nikov Their calculations were based primarily on 
geometric arguments and the common properties of the stress-energy tensor rather than on a field theory. Using the 
methods of quantum field theory the expressions for (ip 2 ) and (T^) of a scalar field in static spherically symmetric 
asymptotically flat spacetimes have been obtained by Anderson, Hiscock, and Samuel [4j]. They assumed that the field 
is massive or massless with an arbitrary coupling £ to the scalar curvature and in a zero temperature quantum state 
or a nonzero temperature thermal state. The result was presented as a sum of two parts, numerical and analytical: 

numeric analytic • (2) 

The analytical part of their expression is conserved. This has a trace equal to the trace anomaly for the conformally 
invariant field. For these reasons they proposed to use (T{f) ana i yt i C directly as an approximation for (T^) ren . An 
analogous result has been obtained by Groves, Anderson, and Carlson [l^ in the case of a massless spin | field in a 
general static spherically symmetric spacetime. 

There are two questions about all of these approximations. First, what are the limits of applicability of these 
approximations? And second, how can one describe quantum states not considered in the above mentioned works? 

The DeWitt-Schwinger expansion is independent of the quantum state. The criterion for the validity of the DeWitt- 
Schwinger approximation is well known: mL 3> 1. The validity of the other approximations discussed above was 
investigated by the authors by means of comparison with earlier known results. 

In this paper, approximate expressions for ( ( p 2 ) ren and (T{f) of a quantized scalar field in static spherically 
symmetric spacetimes are derived. The field is assumed to be both massless or massive with an arbitrary coupling £ 
to the scalar curvature R, and in a zero temperature vacuum state. The expressions for (<p 2 ) ren and (T>f) ren are divided 
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into low- and high-frequency parts. The Anderson-Hiscock-Samuel approach jj] is used for derivation of the high- 
frequency contributions to these quantities. Being an ultraviolet effect, the high-frequency contribution is general in the 
sense that its value does not depend on the quantum state in which (<p 2 ),. e „ and (Tjf) ren are taken. This contribution 
contains all the ultraviolet divergences and can be renormalized. The low-frequency contributions are determined 
by the quantum state and, as an example, such contributions are calculated in asymptotically flat spacctimc in a 
quantum state corresponding to the Minkowski vacuum (in generally accepted terminology this corresponds to the 
choice of the Boulware quantum state). Both parts of (T£) are separately conserved. For a conformally invariant 
field the trace of the high-frequency part of (TJf) is equal to the trace anomaly and that of the low-frequency part 
of {Tj}) ren is equal to zero. The validity of the approximations is discussed. 

In Sec. II the expressions for the Euclidian Green's function and the unrenormalized stress-energy tensor of a scalar 
field with arbitrary mass and curvature coupling in a general static spherically symmetric spacetime are derived. In 
Sec. Ill the WKB approximations for (p> 2 ) and (T^) of a very massive field are discussed. Section IV describes 
the method of calculating the high-frequency contributions to these quantities and the renormalization procedure for 
(ip 2 ) and (T{f). The low-frequency contributions in the case of asymptotically flat spacetimes are derived and the 
renormalized expressions for (ip 2 ) and (T/, 1 ) are displayed in Sec. V. The results are summarized in Sec. VI. In the 
Appendixes some technical results are derived and long expressions are displayed. 

The units h = c = 1 are used throughout the paper. 



II. AN UNRENORMALIZED EXPRESSION FOR (p 2 ) AND (T$) 



In this section the main points of the Anderson-Hiscock-Samuel approach [4j for obtaining unrenormalized expres- 
sions for (p 2 ) and (T») are outlined. 

First of all the metric of the static spherically symmetric spacetime is analytically continued into Euclidean form 



ds 2 = fch 2 + dp 2 + r 2 (d9 2 + sin 2 9 dip 2 ), 



(3) 



where / and r are functions of p, and r is the Euclidean time (r = it, where t is the coordinate corresponding to the 
timelike Killing vector, which always exists in static spacetime). 

The regularization of (p 2 ) and (T^) is achieved by using the method of point splitting. When the points are 
separated one can show that 



(f 2 ) unren— G e(x , x) , 



(4) 



(T^) unre „=(l/2 - W*G E .,av + 9vG E f & ) + (2£ - 1/2)^. 9 ™G B;ct5 - Z(G E f v 

+9^iG E , & p) + 2^(m 2 + £R)G E + t(R£ - 5»R/2)G E - <5> 2 G B /2, 



(5) 



where m is the mass of the scalar field, £ is its coupling to the scalar curvature R, and is the bivector of parallel 
transport of a vector at x to one at x. 

The integral representation for the Euclidean Green's function G E (x,x) of a scalar field in a static spherically 
symmetric spacetime used by Anderson et al. Q is the following: 



G E (x; x) 



— / dwcos[o;(r-f)]^(2Z + l)Pi(cos7) C ui Pui{p<) QW(P>): 

n ^° (=0 



(6) 



where Pi is a Legendre polynomial, cos 7 = cos 9 cos 9 + sin 9 sin 9 cos(p — p) , C^i is a normalization constant, p < and 
p> represent the lesser and greater of p and p, respectively, and the modes p u i (p) and q^i (p) satisfy the equation 
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and the Wronskian condition 
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Above, it is assumed that the field is in a zero temperature vacuum state defined with respect to the timelike Killing 
vector. 



By the change of variables 



qui 



V2r 
1 



i=exp{/V r V^}, 



V2r 2 W 



exp 



one sees that the Wronskian condition ||S} is satisfied by C^i = 1 and the mode equation J7J) gives the 
equation for W(p): 

TT , 2 2 /,„ ^ / x , /' (T4^ 2 )' /(T^ 2 )" 5/(W 2 )' 2 

where 



VWm 2 + (2^UV (r2) " ' /V) ' ^ 



4 / r 2 
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2/ 2 2r 4 fr 



The prime denotes a derivative with respect to p. 

Now one can rewrite expressions (@J,© using expressions J^J,© and then suppose p = PjO 
superficial divergences in the sums over I that appear in this case can be removed as in Refs. |2j, 
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III. WKB APPROXIMATION FOR {<^ 2 ) AND (T^) OF A VERY MASSIVE FIELD 

Obtaining the exact solution of Eq. 110(1 is a very complicated problem. However, it can be solved iteratively if 
there is a small parameter of the considered task. 

Let us evaluate the order of terms in this equation. If the characteristic scale of variation of the gravitational field 
is designated as L: 



L-^maxjll/rl^lnC/r 2 )]'!, | [ln(/r 2 )] " | V2 , |[ln(/r 2 )]' 
then in case r c = 1/m <C £, one can consider the quantity 



1/3 



(21) 



(22) 



as a small parameter of the iterative procedure. In this case as the zeroth-order term of the iterative procedure one 
can choose 



(W 2 ) [Q) = LU 2 + ^l{l + l) + fm 2 



(23) 



-2 , , 2 1 / id 1 ~\\ 1 -p„,2 

The term of the next order [the iteration procedure demands the using of (W 2 )^ instead of W 2 in the right hand 
side of Eq. 1)1 OJl for calculation of the next order term] 



/ , t t , f(W 2 ) {0) ' , f(W 2 ) (0) " 5f(W 2 ) (0) ' 



8 (W 2 ) {0) 4 (W 2 ) {0) 16 (w 2 ) 
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has the order (£ W kb) 2 (W /2 )(o)j i- e -> f° r & H values of w and I this term is much less than the zeroth-order term. For 
derivation it is convenient to choose as the zeroth-order term of the iterative procedure the following: 

(W 2 ) i0) =uj 2 + t(l + ^j + fm 2 , (25) 

because the addition to (VF 2 )( ) some of the addends from (W 2 )( 2 ) does not change the conclusions: 

(W 2 ) {0} » (W 2 ) (2) » (W 2 \ A) » . . . (26) 

or 

(^ 2 )(2)~4kb(^ 2 ) ( o), (^ 2 ) (4) ~4 kb (^ 2 ) (0) , ... (27) 

and 

W 2 = (W 2 ) {0) + {W 2 ) (2) + (W 2 \ 4) + .... (28) 

This approach gives (Ref. 0) the DeWitt-Schwinger expansions of (ip 2 ) and (Tjf) in terms of powers of l/(mL). Note 
that there are the different methods of derivation of the approximations for (ip 2 ) and (T^) in the large mass limit [? 



IV. HIGH-FREQUENCY CONTRIBUTION TO (ip 2 ) AND (Tff) 



In the case r c > L or for a massless field a small parameter does not exist. But we can evaluate the contribution 
to (T^) of the high-frequency modes. For that it is necessary to impose a lower-limit cutoff wq on the integrals over 
W 



in expressions {TBJ-JTSJ {u > uq = way/r 2 / f). As a small parameter of the iterative procedure we can choose 
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This also means [see the definitions of L, Eq. (|21|) . and u, Eq. I|2(J|) ] 

» ^ > 1. 

The zeroth-order solution of Eq. I|1(J|) can be chosen as follows: 

(W 2 ) {0) = 

Then the second and fourth orders are 
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The high-frequency contributions to the quantities B\, B^ B$, and B^ are obtained by substituting the WKB 
expansion of W 2 into expressions (|16|) - (|19[) and imposing a lower-limit cutoff uq on the integrals over u: 
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The mode sums and the integrals in these expressions are of the form 
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where k and n are integers, k > and n > —1. The subtraction terms for the sum over I can be obtained by 
expanding the function being summed in inverse powers of / and truncating at 0(1°). This subtraction corresponds 
to the removing of the superficial divergences in the sums over I discussed above: 
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S"(£,uo) = / du cos(eu) 
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For the other quantities S*(e, Uo) there are no subtraction terms. The details of the calculations of 5*(e, tto) in the 
limit e — > are discussed in Appendix A: 
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The substitution of these expressions into Eqs. an d then into Eqs. fT3" |l -lfT5 |l gives (<p 2 )unren an d 

■" '" r ''' ,. Let us note that the expansions of ((p 2 )unren an d 



the high-frequency contributions to (ip ) unre „ and (T^) 
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in terms of powers of uq correspond to the DeWitt-Schwinger expansions of ((p 



and (TP) 

unren 



terms of powers of ml. 

The renormalizations of (ip 2 ) and {T^) are achieved by subtracting the renormalization counterterms from (ip 2 
and {Tj^) unren and then letting f — > t: 
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C is Euler's constant, R a p is the Ricci tensor, R = R™, a is one-half the square of the distance between the points 
x and x along the shortest geodesic connecting them, is the covariant derivative of a, and the constant m DS is 
equal to the mass m of the field for a massive scalar field. For a massless field m DS is an arbitrary parameter due to 
the infrared cutoff in (i^ 2 )^. A particular choice of the value of m DS corresponds to a finite renormalization of the 
coefficients of terms in the gravitational Lagrangian and must be fixed by experiment or observation. The expression 
for (T^) Dg is displayed in Ref. and is too long to be displayed here. 

For the metric © the calculations of quantities and are similar to those in Ref. p|: 
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The high-frequency contributions to (<f 2 ) and {T^) unren contain all the ultraviolet divergences and can also be 
renormalized. This procedure gives the second-order WKB approximation for { l f 2 ) ren and the fourth-order one for 
components of (Ttf) ren that contain the contributions of high-frequency modes only (w > wq): 
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where (t/? 2 )^, (<P 2 )( 2 v anc ^ tne nontrivial components of (T^) WKB of zeroth and second WKB orders have the form 
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The quantities of the second WKB order (jp 2 ), 2 ) an< ^ ^ ne f° ur th WKB order {Tjf)^ are equivalent to the analytical 
approximations (f 2 ) ana i ytic an d ^Xv) analytic °f Anderson, Hiscock, and Samuel Q (the constants m DS and wq are 
equal to fi and A in their expressions for (<P 2 ) analytic an d (T£) analytic respectively). In the coordinates © the 
expressions for (T>f)^ are given in Appendix B. 

The quantity (T^) WKB is conserved and has the trace 



\ w A 



M 



~2~ 



'9/' 2 3/'(r 2 )' 3/"' 



4/ 2 



2/r 



2/ 



M/(4) 



(69) 



For a conformally invariant field this trace is equal to the trace anomaly. 



V. LOW-FREQUENCY CONTRIBUTION TO (tp 2 ) AND (T£) 

The behavior of low-frequency modes is determined by the boundary conditions and the topological structure of the 
spacetime. If the spacetime is asymptotically flat and the characteristic scale of the gravitational field inhomogencity 
A is much less than the parameter y/J / ujq , 

A « 1, (70) 



V7M) 

the low-frequency contributions to (ip 2 ) and (T^) can be expanded in terms of powers of this small parameter. In 
the following the zeroth term of this expansion will be used for approximation of the low-frequency contributions to 
(ip 2 ) and (T^). This means that we choose the long- wave modes approximately coincident with long- wave modes of 
the Minkowski vacuum (in generally accepted terminology this corresponds to the choice of the Boulware quantum 
state). For these modes [ds 2 = dT 2 + dx 2 + dy 2 + dz 2 = dT 2 + dr 2 + r 2 (d8 2 + sin 2 dip 2 )] 

1 [ Jnj3 exp (iflAT + ip a Ax a ) 



Ge{x, x) — — — — r / dttd A p 



(2tt) 4 J " (tt 2 + pi + p 2 , + p 2 z + m 2 
1 fjnJiAT psm(pAr) 



dVLe 1 ^ 1 / dp 



Air 3 J J Q Ar (tt 2 + p 2 + m 2 ) 

1 f , ,sjat ex P (- AryiF" 
"8^ 2 ' 
1 

: 8^ 



dClt 



iflAT 



Ar 



O(Ar) 



(71) 
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The first addend in the subintegral function must be removed because it gives a superficial divergence similar to that 
discussed above. Then the low- frequency contribution to (<p 2 ) is 



(v? 2 )lfc= lim 



At^O 8tt 2 



dSle inA V^ 2 + m 2 



^o\/W + m 2 +m 2 In 



O + V^o 2 +m 2 



to 



In the case of a massless field or if f2 3> to this expression can be rewritten as 



(72) 



/ 2\ ^0 



8tt 2 16tt 2 16tt 2 
If we take into account f2g = w o/ V 7 / then 

(<£ 2 ) re „=lim [(^Wen ~ (^ 2 )ds] = lim [(f 2 ) 

4«o 2 



In 



4r» 2 



O 



T — T 

R 



HFC 

unren 



(<P 2 )ds] + (^ 2 )lFC = (^ 2 )wKB + (^ 2 )lFC 



TO 2 s r 2 



■f 2 r . f(r 2 )' 

96tt 2 / 2 96tt 2 / 967r 2 /r 2 



The corresponding expressions for (^) are 



4 4 
TO* TO , 

+ —^t In 



LFC 16tt 2 16tt 2 128tt 2 64tt 2 



m 2 



TYl 



L 2 u^ 



o\ — 



(73) 



(74) 



(75) 



<^> LFC =(^> LFC = {T$\ 



+ 



+ 



3m 



+ 



■In 



m 2 



TO 



V/lfc 487r 2 lg7r 2 1287r 2 647r 2 

The analytical approximation for (T^) ren in the case of asymptotically flat spacetimes is 



o\ — 
1 



z> / ren 

lim f(T' 1 ) - (T p ) 1 = lim |"(T^ HFC - (T M ) 1 + (T M ) 

T — >T T — >T 

\ I//WKB 1 \ 1//LFC' 



(76) 



(77) 



( T t) 

\ z I r 



1\ ( 1 5/' 2 f (r 2 )' (r 2 )' f" (r 2 f 
? 67 I 2r 2 4/ 2 /r 2 8r 4 / 2r 2 



4 4 
TO* TO* , 

+ — — In 



128tt 2 64tt 2 



4« 2 



(2?) 



(78) 



( T p) = 

\ r / ren 



u 



P I ren /^^2y,2 



? 6/ 1 2r 2 4/ 2 4./Y 2 8r 4 



+ 



3to 4 



TO* 

+ 77-^hi 



128tt 2 64tt 2 



4u 2 



+ ( r p P )(4)+ l „2 



1 



U 2 L 4 



\rL 4 



(79) 



PI ren 4^2 1$ g 



3to 4 

+ 



(r 2 )" /" , 5/' 2 (r 2 )' 2 /'(r 2 )' S 



128tt 2 64tt 2 



In 



4 U 2 



4 r 2 



4/ 8/ 2 
O' 1 



8r 4 8/r 2 



ri 4 



(80) 



This tensor is conserved, 



ren; u 



O 



U §L« 



(81) 
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and has the trace 



3/V/_ 3T S 
2/ 




16tt 2 



16tt 2 



■In 



An 2 



For a conformally invariant scalar field this trace is also equal to the trace anomaly. 



VI. CONCLUSION 



(82) 



In this paper, analytical approximations for (<p 2 ) r en and (T*f) ren of quantized scalar fields in static spherically 
symmetric asymptotically flat spacetimes have been obtained. The field is assumed to be in a zero temperature vacuum 
state, with mass m<l/L \L(p) is the characteristic scale of variation of the gravitational field at the considered point] 
and with an arbitrary coupling £ to the scalar curvature. 

The necessary conditions for the validity of the analytical approximations H74|) . l|78[) - H80|) are 



A « « L{p) f 



Wq 



(83) 



where A is the characteristic scale of the gravitational field inhomogeneity in asymptotically flat spacetime and Wq is 
the constant of the WKB expansion. 

If we consider a spherical body with radius r$ > r g {r g is the gravitational radius of this body, A ~ ro) then the 
metric spacetime outside this body is 



ds=-\l- 



r{p) 



dt 2 + dp 2 + r(p) 2 (d6 2 + sin 2 9 dip 2 ), 



(84) 



where r{p) is the inverse function to the function 



p(r)=p + r g 



■In 



r„ \r„ J \ r„ \ r„ 









- 1 - V r o/ r g) 


(Vr/r g - 1 - 


- V r / r 9 ) 




- 1 + V r o/r g ) 



When r(p) ^> r 

/(/0~1, L(p)~r(j>), 
and the conditions (|83|l can be satisfied by the choice of wq, 

r < w < r(p). 



(85) 



(86) 



(87) 



This means that, far from the body where r(p) ^> ro, the approximations (|74|l . (|78|l - (|80|l are valid. 

The presence in expression s II74II.(|78I) - (I80|I the arbitrary parameter uq = wo r /v / 7 ^ s a generic feature of local 
approximation schemes 0, Il2l Il5l Il8l Il3| . For a conformally invariant field this parameter can be absorbed into the 
definition of the constant m DS . 

When the massless conformally coupled scalar field is in a zero temperature vacuum state and propagating in static 
spherically symmetric asymptotically flat spacetimes the expressions H74|) . l|78|l - H80() are equivalent to the approxima- 
tions of Page, Brown, and Ottewill 0, IE E| , Zannias , and Frolov and Zel'nikov ^l| (for the particular choice of 
the arbitrary parameters in their expressions), and the analytical approximation of Anderson, Hiscock, and Samuel 
0. Let us note that in this case the low-frequency contributions to ((p 2 ) re n and {Tjf) are equivalent to the low- 
frequency contributions which are given by the procedure of |4| for modes with w < wq. This means that using such 
calculations in pol l2ll l22l | is correct in the case of a conformally invariant field. Let us note also that the more exact 
procedure [23| gives different result. 
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Appendix A 

To calculate the quantities S^(e,uo) it is necessary to compute the various sums over I. We start from the sum in 
expression (l3H : 



S(u,n) = 



+ 1) 



2n+l 



i 2 + (Z + l/2)^]( 2 ™- 1 )/ 2 



For calculation of this sum we can use the Abel-Plana method llfll 124. |2 



l + - 2 ) +(2n-l) y 



q+lo ° F(z) 



n > 0. 



ou ,.00 W(z\ rq+ico 

E m + 1/2) = / F{ X )dx + / Y-^Ldz / - 



dz, 



where —1/2 < q < 1/2, f(z) is a holomorphic function for Rez > q, f(z) satisfies the condition 

\F(x + iy)\ < e(x)e aM , < a < 2vr, 
and e(:r) is an arbitrary function with asymptotic behavior 

e(x) — ► for x — ► +oo. 
Using this formula we can calculate the sum IjAlfl : 



S(u, n)— lim 



r 2n+l 



g+zoo 



(u 2 +a; 2)(2„-l)/2 
r 2ri+l 



. 2 + (2n-l)^- 



( u 2 + z 2)(2n-l)/ 



-_-^ + (2n-l)- 



v 2n+l 



_^ + (2n-l) 



( u 2 + z 2)(2n-l)/ 



(2,-l)(2n)!!u3 m 
3 (2n- 1)!! v ' s^+o 



(terms of this integral that 
diverge in the limit S — > +0 

(2n- 1) u 3 j (2n)!! / rf 
(2ra-l)!! j - 3 Vwf^ 



U+<5 



(1 + e i2,rz ) 

(1 + e - i2nz ) 
x 2n+1 dx 



o (^-^^"-^(l + e 2 ^) 



1 " 
/ 1 + e 27 ™?> 

The sum in expression l|40[) can be calculated by the differentiation of S(u,n): 



(n > 0). 



^\[u 2 + (l + l/2) 2 }( 2n + 1 )/ 2 
(2n)!! u 2 ( d 



-1 = 



-1 



(2n — 1) \udu 



S(u, n) 



,3 M 



„ T g/vj - y 2 ^ 

7 1 + e 2 ™?' 



(Al) 

(A2) 

(A3) 
(A4) 



(A5) 



(n>0). (A6) 



(2n-l)!! (2n - 1)!! \udu J 1 V M du, 
The other sums in the expression for S*(e,u ) [Eq. (J3SJ] can also be calculated by the differentiation of S(u, n): 
{l + h) 2n+1 _ (-l) fc (2n-l)!! ( d\ k+2 



E 



£ kt 2 + + l/2) 2 ](2«+2fe+3)/2 (2n - l)(2n + 2fe + 1)!! V^w J 



S(u, n) 
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(2n)!! (2k- 1)!! 

M 2fc+1 

fc+2 I 



(2n + 2A: - 
2(-l) fc 



1)!! 



f — 

(2n + 2/c + l)!! 



f— Y 

\udu J 



,2)1 



1 ?/\A - y 2 rfy 
1 + e 27Tu y 



n > 0, 
fe > 



(A7) 



The integral over y in all these expressions can be calculated approximately. In this integral the integrand decreases 
exponentially for y > 1/(2ttu). The parameter u in all expressions l|38(l satisfies the condition u > uq 3> 1. Therefore 
the main contribution to this integral gives a small value of y (y <C 1) and the square root in the integrand can be 
expanded. Then the up limit can be gone to infinity: 



iV l - y 2 dy 

I _|_ f,2nuy 



ydy 



V 2 1!! 4 3!! 6 



i 



7 



31 



48u 2 3840w 4 129024m 6 



(A8) 



Substituting this expression into Eqs. ljA5(l - l|A7(l and integrating in Eq. I|38|) we obtain the resulting expressions 



Appendix B 



T 



( 4 ) 460807r 2 r 8 /' 



|360r 8 / 4 ™ 4 + 360r 8 / 2 /' 2 - 480T 6 / 3 /'(r 2 )' - 480r 8 / 3 /'' 



+32r 4 / 4 - 32(r 2 )""r 6 / 4 - 12/" 2 r 8 / 2 + 56(r 2 )"V/ 4 + 16/""r 8 / 3 + 7/'V 



-32/" (r 2 )*V b / 3 - 16/V)'V/ 3 + 40/' 2 (^)* V e '/ a - 16ff"r»f + 16/"(r 2 )*V/ 3 
+4/' 2 /"r 8 / - 12/' 3 (r 2 )V 6 / - 14/' 2 (r 2 )' 2 r 4 / 2 - 48/'(r 2 )"V 6 / 3 + 48(r 2 )'(r 2 )"V 4 / 4 
-112(r 2 )' 2 (r 2 )"r 2 / 4 + 32/'(r 2 )"(r 2 )V 4 / 3 + 40(r 2 )' 4 / 4 + 32(r 2 )'/"'r 6 / 3 



5760r 4 / 4 - 11520(r 2 )"r 4 / 4 + 2880(r 2 )'V/ 4 + 8640/ ,2 r 6 / 2 



11520/V/ 3 



+5760(r 2 )" 2 r 4 / 4 + 4320/" 2 r 8 / 2 + 7560/' 4 r 8 + 360(r 2 )' 4 / 4 + 11520/"(r 2 )"r 6 / 3 
+2880/'(r 2 )' 3 r 2 / 3 - 2880/ ,2 (r 2 )"r 6 / 2 + 5760/'/"V 8 / 2 - 2880/"(r 2 )' 2 r 4 / 3 
-17280/ ,2 /"r 8 / - 11520/' 3 (r 2 )V 6 / + 720/ ,2 (r 2 )'V/ 2 + 11520/'(r 2 )"V 6 / 3 
-2880(r 2 )' 2 (r 2 )"r 2 / 4 + 14400/'/"(r 2 )V 6 / 2 - 5760/'(r 2 )"(r 2 )V 4 / 3 l (£ - 1/6) 2 



2880r 6 / 4 - 5760r 8 / 3 /" - 5760r 6 / 3 /'(r 



4320r 8 / 2 /' 2 + 720r 4 / 4 (r 2 



s'2 



2 r4 



-2880r b / 4 (r 2 ) J m 2 - 960/"r b / 3 + 960(r 2 ) V 4 / 4 + 720/'V/ 2 - 960(r 2 ) r 2 f 
-960(r 2 )' 4 / 4 + 1920(r 2 )'/"V 6 / 3 + 960(r 2 )"V/ 4 - 2160/" 2 r 8 / 2 - 1440(r 2 )" 2 r 4 / 4 



+960/""r 8 / 3 + 2880/"(r 2 )"r 6 / 3 + 1200/'(r 2 )' 3 r 2 / 3 - 2160/' 2 (r 2 )"r 6 / 2 



\r 2 fr 2 f 

2880/'/"V 8 / 2 - 1200/"(r 2 )'V/ 3 + 7920/ /2 /"r 8 / + 2640/ /3 (r 2 )V 6 / 
900/' 2 (r 2 )'V/ 2 + 1920/'(r 2 )"V 6 / 3 + 2880(r 2 )' 2 (r 2 )"r 2 / 4 - 3780/' 4 r 8 
1440(r 2 )'(r 2 )"'r 4 / 4 - 2400/'(r 2 )"(r 2 )V 4 / 3 - 4320/7"(r 2 )V 6 / 2 l (£ - 1/6) 



-720m 4 r 8 / 4 - 16r 4 / 4 - 20(r 2 )' 4 / 4 + 56(r 2 )' 2 (r 2 )"r 2 / 4 + 8/'(r 2 )"(r 2 )V 4 / 3 

-28(r 2 )" 2 r 4 / 4 - 16/""r 8 / 3 + 52/'/"(r 2 )V 6 / 2 + 16(r 2 )"V/ 4 + 36/" 2 r 8 / 2 
+49/' 4 r 8 - 4/'(r 2 )' 3 r 2 / 3 - 4/' 2 (r 2 )"r 6 / 2 + 48/'/"V 8 / 2 + 4/"(r 2 )'V/ 3 
-116/' 2 /'V 8 / - 22/' 3 (r 2 )V b / - 3/' 2 (r 2 ) V/ 2 + 8/'(r 2 )'V/ 3 - 32(r 2 )'/"V 6 / 3 
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-24(r 2 )'(r 2 )'"r 4 / 4 - 8f"(r 2 )"r 6 f + (-2880r 4 / 4 + 4320(r 2 )V/ 4 + 5580(r 2 )' 4 / 4 
+5760/'(r 2 )"(r 2 )V 4 / 3 - 15840(r 2 )' 2 (r 2 )"r 2 / 4 + 6480/" 2 r 8 / 2 + 5760(r 2 )"V/ 4 



-2880/"V/ 3 - 5760(r 2 )""r 6 / 4 - 2880/' (r 2 )V/ 3 + 5760/' 2 (r 2 )V/ 2 
+15840/'/" (r 2 )V 6 / 2 + 8640/7'V/ 2 + 2880/"(r 2 )' 2 r 4 / 3 - 20880/' 2 /"r 8 / 
-9360/' 3 (r 2 )V 6 / - 2160/' 2 (r 2 )'V/ 2 - 2880/'(r 2 )"V 6 / 3 - 5760(r 2 )'/"V 6 / 3 
+8640(r 2 )'(r 2 )"V 4 / 4 - 5760/"(r 2 )"r 6 / 3 + 8820/ ,4 r 8 ) (£ - 1/6) 2 



+ (2880r 6 / 4 (r 2 )" - 720r 4 / 4 (r 2 )' 2 - 2880r 6 / 4 ) (f - 1/6) m 2 



In 



(Bl) 



{-1080m 4 r 8 / 4 - 120/' 2 m 2 r 8 / 2 + /' 4 r 8 + 32/'(r 2 )"(r 2 )V 4 / 3 



J 2 



^ 3 



p; ( 4 ) 460807r 2 r 8 / 4 

-24/' 2 (r 2 )V/ 2 + 12/'V)'V/ 2 + 8f'f"'r s f 4f" 2 r s f + 16/"(r 2 ) V/ 3 
-8/' 3 (r 2 )V 6 / + 32/7" (r 2 )V 6 / 2 + [5760/'(r 2 )"(r 2 )V 4 / 3 - 1440/'(r 2 )' 3 r 2 / 3 



-4/'7'V 8 / - 24/"(r 2 f r 4 / 3 - 16/'(r 2 )"V 6 / 3 - 16(r 2 )'/"V 6 / 3 - 8/'(r 2 f r 2 / 3 

'I' 



,2 



-5760/V) V 4 / 3 - 2880/' V/ 2 + 2880/7'V) > / + 2160/'> 2 ) r 4 / 



-720/'V 



1440/ ,2 /"r 8 / + 2880/' 2 (r 2 )"r 6 / 2 l (£ - 1/6) 2 



+ 



240/' 2 r 6 / 2 + (720(r 2 )'V/ 4 - 1440/'(r 2 )V 6 / 3 - 1440/' 2 r 8 / 2 

-2880r 6 / 4 ) m 2 - 1260/' 2 (r 2 )'V/ 2 + 540/'V + 480(r 2 )'/"V 6 / 3 - 720/' 2 /"r 8 / 
+480/7"'r 8 / 2 - 480/"(r 2 )"r 6 / 3 - 720/'(r 2 )"(r 2 )V 4 / 3 - 1200/'/"(r 2 )V 6 / 2 



+1200/ ,2 (r 2 )"r 6 / 2 + 480/'(r 2 )'V/ 3 + 720/'°(r 2 )V t) / + 720/"(r 2 )'V/ 3 
+480/'(r 2 )V 4 / 3 - 240/" 2 r 8 / 2 ] (£ - 1/6) + [-720m 4 r 8 / 4 + 12/"(r 2 )' 2 r 4 / 3 

-16r 4 / 4 - 4(r 2 )' 4 / 4 + 8/' 2 (r 2 )"r 6 / 2 + 12/ ,2 /"r 8 / + 8(r 2 )'/"V 6 / 3 - 8/7"'r 8 / 2 
+8/'(r 2 )'V/ 3 - 24/'/"(r 2 )V 6 / 2 - 3/' 2 (r 2 )' 2 r 4 / 2 + 4/'(r 2 )' 3 r 2 / 3 
-16/'(r 2 )"(r 2 )V 4 / 3 + 8(r 2 )' 2 (r 2 )"r 2 / 4 - 8(r 2 )'(r 2 )"V 4 / 4 + 10/' 3 (r 2 )V 6 / 
-7/' 4 r 8 - 8/"(r 2 )V/ 3 + 4/" 2 r 8 / 2 + 4(r 2 )"V/ 4 + (-2880r 4 / 4 

-4320(r 2 )' 2 r 2 / 4 - 2880/' 2 (r 2 )"r 6 / 2 - 2340(r 2 )' 4 / 4 - 1440/7"'r 8 / 2 
+2880/"(r 2 )"r 6 / 3 + 4320/'/"(r 2 )V 6 / 2 - 1440/' 3 (r 2 )'r 6 f + 4320f 2 (r 2 fr\f 2 
+720 f(r 2 fr 2 .f ' 3 - 2880/'(r 2 )"V 6 / 3 + 5760(r 2 )' 2 (r 2 )"r 2 / 4 - 5760(r 2 )'(r 2 )"'r 4 / 4 



+lU0f'(r 2 )"(r z )'r\f - 1260/'V - 2880(^)'/'"r e / 3 + 2160/'7"r 8 / 
-4320/"(r 2 )' 2 r 4 / 3 + 720/" 2 r 8 / 2 + 2880(r 2 )"V/ 4 ) (£ - 1/6) 2 + (-2880r 6 / 



„2\' till Jo j-3 



c' 2 f"„8 . 



6 j-4 



+720(r 2 )' 2 r 4 / 4 + 1440/'(r 2 )V 6 / 3 ) (£ - 1/6) m 2 



In 



4u 2 



m 2 DS r 2 



(B2) 



W) (4 ) = ( T < 



1 



(4) 460807r 2 r 8 / 4 



{-1080r 8 / 4 m 4 + 120r 8 / 2 / ,2 m 2 - 8r 4 / 3 /'(r 2 )"(r 2 )' 



+ 16r 6 / 2 /'/"(r 2 )' + 17r 8 /' 4 - 16r 8 / 3 /"" - 8r 6 //' 3 (r 2 )' + 16r 6 / 3 /'(r 2 )"' 
+ 16r 4 / 3 /"(r 2 )' 2 - 52r 8 //' 2 /" + 24r 8 / 2 /'/"' - 4r 4 / 2 /' 2 (r 2 )' 2 
+8r 6 / 2 /' 2 (r 2 )" + 28r 8 / 2 /" 2 - 24r 6 / 3 (r 2 )7"' - 16r 6 / 3 /"(r 2 )" 
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-5760r 6 /7" ~ 1440r 4 / 3 /"(^ 2 )' 2 + 2880r 6 / 2 /' 2 - 14400r 4 f f \r 2 )" \r 2 )' 
+ 11520r 6 / 3 /'(^)"' + 2880r 8 / 2 /" 2 + 5760r 6 f f" \r 2 )" - 7920r 6 / f (r 2 )' 
-3600r 4 /7'V)' 2 - 14400r 8 //' 2 /" + 8640r 4 / 3 /'(^ 2 )' + 2880r 6 / 2 .f' V)" 
+5040r 2 / 3 /'(r 2 )' 3 + 6480r 8 /' 4 + 5760r 8 / 2 /'/"' + 10080r 6 / 2 /7"(^ 2 )'l (£ - 1/6) 2 

-240r 6 / 2 /' 2 + 480r 6 / 3 /" + (-720r 4 / 4 (r 2 )' 2 - 1440r 8 / 3 /" + 720r 8 / 2 /' 2 
+1440r 6 / 4 (r 2 )" - 720r 6 / 3 /'(r 2 )') m 2 - 1440r 8 / 2 /" 2 - 720r 4 / 3 /' (r 2 )"(r 2 )' 

-1920r 8 / 2 /7"' + 5400r 8 //' 2 /" - 2580r 8 /' 4 + 720r 6 / 3 /"(r- 2 )" + 480r 8 / 3 /"" 
-2280r 6 / 2 /'/"(r 2 )' + 720r 6 / 3 (r 2 )' - 720r 4 / 3 /"(r 2 )' 2 + 1440r 6 //' 3 (r 2 )' 
+480r 2 / 3 /'(r 2 )' 3 + 660r 4 / 2 /' 2 (r 2 )' 2 - 960r 6 f 2 f (r 2 )"] (£ - 1/6) 

-720r 8 / 4 m 4 + 16r 4 / 4 - 20r 8 ff'f" + 12/ 4 (r 2 )' 4 - 14r 6 / 2 /'/"(r 2 )' 

+12r 8 /V)7"' + 8r 6 / 3 / V)" + ^ f f (r 2 )' 2r e f 2 f' 2 (r 2 f - 8r 6 /W)'" 
+ 12r 4 / 4 (r 2 )" 2 - 20r 8 / 2 /" 2 + 4r 4 / 3 /'(r 2 )"(r 2 )' - 8r 4 / 3 /"(r 2 )' 2 + 3r 4 / 2 /' 2 (r 2 )' 2 
+ 16r 4 / 4 (r 2 )'(r 2 )"' + 52r 8 //' 2 /" + 8r 8 / 3 /"" - 8r 6 /V 2 )"" - 32r 2 / 4 (r 2 )' 2 (r 2 )" 
+ (2880r 4 / 4 - 21r 8 /' 4 - 8640r 4 / 4 (r 2 )" - 4320r 4 / 3 /'(r 2 )' + 8640r 2 / 4 (r 2 )' 2 

+9360r 6 / 2 /7"(r 2 )' + 7200r 8 / 2 /7"' - 3240r 4 / 2 / ,2 (r 2 )' 2 + 5040r 6 / 2 /' 2 (^ 2 )" 
+7020/ 4 (r 2 )' 4 - 5400r 2 / 3 /'(r 2 )" 3 - 20880r 2 / 4 (r 2 )'V)" + 11520r 4 / 4 (r 2 )'(r 2 )"' 
-16560r 8 //' 2 /" - 5400r 6 //' 3 (r 2 )' + 8640r 4 / 4 (r 2 )" 2 + 11520r 4 / 3 /'(r 2 )"(r 2 )' 
-5760r7 3 /V)"' - 7200r 6 / 3 /"(r 2 )" - 4320r 6 / 3 (r 2 )7"' - 2880r 8 / 3 /"" 
+5040r 8 / 2 /" 2 + 5040r 4 / 3 /"(r 2 )' 2 - 5760r 6 / 4 (r 2 )"" + 7020r 8 /' 4 ) - 1/6) 2 

+ (-720r 4 / 4 (r 2 )' 2 + 720r 6 / 3 /V)' + 1440r 6 / 4 (r 2 )" + 1440r 8 / 3 /" 



-720r 8 / 2 /' 2 ) (C-l/6)m 2 



In 



4m 2 



o 



(B3) 
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